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Abstract. We present the first comprehensive measurement of D°,D + ,D* + and 
^ • their charge conjugate states at mid-rapidity in d+Au collisions at % /s JVJV =200 GeV 

using the STAR TPC. The directly measured open charm multiplicity distribution 
£f) \ covers a broad transverse momentum region of 0< pr < H GeV/c. The measured 

dN/dy at mid-rapidity for D° is 0.0265 ± 0.0036(siai.) ± 0.0071(syst.) and the 
measured D* + / D° and D + /D° ratios are approximately equal with a magnitude of 
' 0A0±0.09(stat.)±0.13(syst.). The total cc cross section per nucleon-nucleon collision 

■ extracted from this study is 1.18±0.21(s£at.)±0.39(sysi.) mb. The direct measurement 
(T^) [ of open charm production is consistent with STAR single electron data. This cross 

section is higher than expectations from PYTHIA and other pQCD calculations. The 
measured px distribution is harder than the pQCD prediction using the Peterson 



o 

' fragmentation function. 

=5 



PACS numbers: 25.75.Dw, 13.25.Ft, 14.40.Lb 



1. Introduction 



X ■ Study of open charm production provides valuable test of perturbative QCD (pQCD) 
predictions for heavy quark production [1]. There are significant uncertainties 
in the leading order (LO) and next-to-leading order (NLO) pQCD calculations 
of charm quark production depending on choices of the charm quark mass, the 
renormalization/factorization scale and the parton distribution function (PDF) [2]. 
Therefore, the open charm measurement provides a constraint on the parameters used in 
the pQCD calculations. The energy dependence of the cc production cross section is an 
important issue. Previously, most measurements of open charm production in hadronic 
interactions were from fixed-target experiments at lower energy (^/s <63 GeV). Only 
recently, CDF reported new results on open charm production in pp collisions in the 
high pt region from the collider experiment at y^s=1.96 TeV [3]. The open charm cross 
section has been studied in Au+Au collisions at y/s NN =130 GeV indirectly through 
measurements of single electrons from the open charm semi-leptonic decay [11]. 

The charm quark is generally believed to be produced primarily through gluon 
fusion (gg — > cc) in parton-parton hard scattering due to its heavy mass, which makes 
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the charm quark an important probe to the properties of matter formed in the early 
stage of relativistic heavy ion collisions. Theoretical calculations predicted that the 
open charm cross section in AA collisions would be significantly enhanced with respect 
to purely hadronic production [5] as a consequence of quark-gluon plasma (QGP) 
formation. Some studies have already attributed the increased yield of dileptons in the 
intermediate mass region reported by NA50 to open charm enhancement [6]. On the 
other hand, the open charm momentum spectrum may also be modified in AA collisions 
due to final state interactions such as parton energy loss in medium. However, it is 
predicted that energy loss of a charm quark should be considerably smaller compared 
to a light quark as a consequence of the dead-cone effect, leading to an enhanced D/n 
ratio as a function of px [7]. 

It is found that the relative fraction of the charm quark hadronization into different 
open charm species is nearly the same in elementary particle collisions ('jp, e + e~, 
ep etc) [8]. Could the charm quark hadronization be modified in AA collisions? A 
calculation based on the Statistical Coalescence Model [18] does show that ratios of 
open charm species in AA collisions would be different from those in the elementary 
collisions. Furthermore, at RHIC, the predicted J ftp suppression has to be investigated 
with respect to the total open charm cross section since the Drell-Yen cross section, 
which is the reference of the J/ip suppression at SPS, is small in comparison with the 
open charm cross section at RHIC. In order to investigate the outstanding physics issues 
mentioned above, it is crucial to study open charm production in a less complicated 
system, where the nuclear effects, like QGP formation, parton energy loss, are not 
expected. The d+Au experiment at RHIC provides such an essential reference. 

2. Data Analysis 

The main detector involved in this data analysis is the STAR Time Projection Chamber 
(TPC) in a 0.5 T solenoidal magnetic field. A Zero Degree Calorimeter (ZDC) in the Au 
beam direction was used for the minimum bias trigger by requiring at least one beam- 
rapidity neutron in this ZDC. This trigger accepts 95±3% of the ci+Au hadronic cross 
section. For a detail description of the STAR detector, see [10]. In this study, about 
15 million d+Au minbias events with primary vertex position in beam direction within 
±75 cm around the TPC center were analyzed. The primary vertex reconstruction 
efficiency is 93± 1% of triggered minimum bias events. The correlation between the 
measured momentum and the ionization energy loss (dE/dx) of charged particles in 
the TPC gas were used for particle identification. The measured (dE/dx) can be well 
described by the Bethe-Bloch function smeared with a resolution of width a. 

The open charm are reconstructed through their hadronic decay channels (and 
their charge conjugates): D° -> K~it + (BR=3.8%), D° -> K-n + p°(BR = 6.2%) -> 
K-Ti+Ti+Ti-, D+ -> ^-7r+7r+(BR=9.1%) and D*+ -> D°n+(BR = 68%) -> K-tt + k+. 
Event mixing technique [12] was used to reconstruct open charm mass peaks. The 
measured charged multiplicity distribution in the region of \r)\ < 0.5 was divided into 
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40 bins and the primary vertex position distribution was divided into 15 primary vertex 
bins. Each event was mixed with 5 other events within the same multiplicity bin and 
primary vertex bin. Tracks within 2 a of the pion Bethe-Bloch curve and 3 a of the 
kaon Bethe-Bloch curve were selected. In the momentum region of pt > 0.6 GeV/c 
where the kaon and pion dE/dx bands merge, the nominal kaon and pion masses were 
assumed in turn for each track. It was also required that the closet approach distance 
of a track to the primary vertex, dca, be smaller than 1.5 cm with exception of the soft 
pions, ir s , from D* decays where dca < 3 cm was used. 

In the following subsections, We will discuss the analysis method for each decay 
channel. 

2.1. D° -> K--K+ (+c.c.) 

For this analysis only tracks within \r}\ < 1 were accepted. The transverse momentum 
and total momentum cut of a track were between 0.2 to 10 GeV/c and 0.3 to 10 GeV/c, 
respectively. The resulting invariant mass distribution after event-mixing background 
subtraction was shown in Fig. (la), where a clear D° signal was seen alone with a residual 
background distribution. We fit the histogram with a Gaussian function for signal plus 
a linear background. From the fit, the mass and width (a) of the D° signal was found 
to be 1.863± 0.003 GeV/c 2 and 13.8±2.8 MeV/c 2 , respectively. 

2.2. D° -> K-7T + p° (+c.c.) 

An independent analysis for D° — > K~7r + p° was performed where only tracks within 
1 77 1 < 1.5 were accepted. The transverse momentum and total momentum cut of a track 
were set between 0.3 to 10 GeV/c. The analysis is similar to that for D° — > K~ir + 
except that additional n + and 7r~ candidate tracks, if their invariant mass satisfied 
0.62 GeV/c < M(7r + 7r~) < 0.86 GeV/c, were combined with the selected K~ and 7r+ 
candidate tracks to form a D° candidate with invariant mass M(K~7r + 7r + 7r _ ). The 
mass plot after event mixing background subtraction is shown in Fig. (lb). We fit the 
histogram with a Gaussian function for signal plus a linear background. From the fit, 
the mass and width (a) of the D° signal were found to be 1.852± 0.005 GeV/c 2 and 
13.6±4 MeV/c 2 , respectively. 

2.3. D* + -> D°ir+ (+c.c) 

The track cuts for kaon and pion candidates from D° decay in the D* + analysis were the 
same as those used in the D° — > K~n + p° analysis. However, special treatments were 
needed for the soft pion daughter, 7r+, which has an average momentum of about 50 MeV. 
For the soft pion candidates it was required that the transverse momentum and total 
momentum cut be set between 0.1 to 1.0 GeV/c and the ratio of the D° to n s momentum, 
p(D°)/p(n s ) > 9.0. First, the invariant mass of kaon and pion candidate tracks, 
M(Ktt), was calculated and the D° candidate which satisfied 1.82< M(Kn) <1.90 
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Figure 1. The D° signal reconstructed 
from D° -> (+c.c.) (a) and 

the -D signal reconstructed from D° — > 

^-7T + p° (+C.C.) (b). 



Figure 2. The £)* signal reconstructed 
from D* — ► i^~7r + 7r+ (+c.c.) (a) and 
the _D + signal reconstructed from D + — > 

X-7T + 7T+ (+C.C.)(&). 



GeV/c 2 were kept. A soft pion candidate, with a charge opposite to that of the 
kaon candidate, was then combined with the D° candidate to form a D* candidate 
with invariant mass M(Kirir s ). Fig. (2a) shows the distribution of the mass difference, 
AM = M(Kir7[ s ) — M(Ktt), for D* ± candidates after all cuts. A clear signal was seen 
around the nominal value of M{D* ± ) — M(D°). We fit the mass distribution with a 
sum of a Gaussian function describing the signal and a linear background. From the 
fit, the mass difference and width (<x) of the D* + signal were found to be 146.37±0.15 
MeV/c 2 and 0.57±0.16 MeV/c 2 . 

2.4. D+ -> K-n + n + (+c.c) 

The track cuts used for the D ± analysis were the same as those for the D* analysis with 
the exception that no special treatment for the pion candidate was needed. Fig. (2b) 
shows the M(Kmr) distribution for the D ± candidates after all the cuts. A clear signal 
was seen at the nominal value of M(D ± ). The mass distribution was fit to a sum of a 
Gaussian function describing the signal and a linear background. From the fit, the mass 
and width (a) of the D* + signal were found to be 1.861± 0.008 GeV/c 2 and 16.9±4.6 
MeV/c 2 . 

3. Results and Discussions 

Fig. (3) shows the invariant multiplicity distributions of the measured open charms as 
a function of pr after the reconstruction efficiency correction, which is about 40%-60% 
(increasing with p T ) except for D* at low p T (eg. 6% at p T =1.5 GeV/c due to the low 
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Figure 3. Measured invariant multiplicity 
distributions for D°, D* + and D + and a 
power-law fit to the data points. 



Figure 4. Consistency between directly 
measured open charm and the open charm 
measurements through single electrons. 



reconstruction efficiency of the soft pions). The spectrum was corrected for vertex and 
trigger efficiency. For D° — > K + tt~, a correction was also done through a Monte Carlo 
simulation for the correlation around the D° mass due to the misidentification of the 
kaon and pion from the D° decay. Errors of the data points are statistical only and the 
data have been divided by 2 assuming a(D) = c{D). Overlap between D° and D* at 
low pt and between the D° and D + at high p T allows us to compare the production 
cross sections of different open charm species. It was found that the inclusive cross 
sections of D* and D + are approximately equal. This observation is consistent with the 
measurements in [3, 8]. 

With the ratio of D* /D°(=D + /D°) as a free parameter we fit the all data 
points by a power-law function, A(l + Pt/Po)~ u , where A,p Q , and n are parameters. 
From the fit, we obtained: dn/dy(D°) = 0.0265 ± 0.0036(stat.) ± 0.0071 (syst.) and 
(p T ) = 1.32 ± 0.08(stat.) ± 0.16(syst.) GeV/c. It was also found that D*/D° ta 
D + 1 D° = 0.40 ± 0.09(stat) ± 0.13(sys). The cc cross section can then be calculated by 
a c5 = 1.24(<tdo +(Td+)i where the factor 1.24 includes the contribution from unmeasured 
D s and charm baryons (A c etc.) [8]. Assuming number of binary collision scaling, one 
could obtain <r D o for nucleon-nucleon collisions by a D o = 4.6 x dn/dy(D°) x cr ine / j 'Nu n , 
where <Ji ne i — 42 mb and Njy in = 7.5 [13] are the inelastic cross section for a nucleon- 
nucleon collision and number of binary collisions in d+A collisions at v /s JVJV =200 GeV, 
respectively. The factor 4.6 accounts for the ratio of the D° multiplicity in the full 
rapidity region to dn/dy(D°) at mid-rapidity and was obtained from a PYHTIA [14] 
calculation with parameters tuned to fit our measured spectrum [15]. With the measured 
dn/dy(D°) and D + /D° ratio, the a c5 for nucleon-nucleon collisions is found to be 
1.18 ± 0.21(stat) ± 0.39(sys) mb. 
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We also studied open charm production through measurements of single electrons 
from the open charm semileptonic decay. Fig. (4) shows the electron spectra for d+Au 
and p + p collisions at v /s^7=200 GeV, after background subtraction in comparison 
with the expected electron spectrum based on the decay of the STAR measured open 
charm spectrum (the bands). The background electrons are mainly from 7-conversion 
in the detector and from the ir° and rj Dalitz decay. The single electrons were measured 
by both Time of Flight (TOF) (solid symbols) and TPC dE/dx (open symbols). The 
width of the bands reflects the uncertainties in the calculation. For details about these 
measurements, see [16]. The open charm measurements through direct constructions of 
open charm signals and through the single electron measurement are consistent within 
STAR. 




Figure 5. Background subtracted single 
electron spectra in p + p and <i+Au 
collisions measured by STAR BEMC 
compared with a PYTHIA calculation. 



Figure 6. The STAR measured D*/D° 
ratio in d+Au collisions compared with 
other experimental data and with theoret- 
ical predictions. 



The single electrons were also measured using the STAR Barrel Electromagnetic 
Calorimeter (EMC) up to 7.0 GeV/c. Fig. (5) gives the background subtracted electron 
spectra for d+Au (circles) and p + p (triangles) collisions together with the calculations 
of PYTHIA with the modified parameters [15]. The PYTHIA curves shows that in 
the region of p T > 4.0 GeV/c the contribution of B meson decays to electrons starts to 
overtake that of D meson decays, which is the clear indication at RHIC that the bottom 
quark contributes significantly to the measured high p? electrons. For details of this 
analysis, see [17]. 

In Fig. (6), the STAR measured D*/D° is compared with the previous 
measurements [3, 8] and with theoretical predictions [18]. Within the experimental 
uncertainties the ratio shows little dependence on collision system and energy. 

The measured open charm spectrum is compared with theoretical predictions from 
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Figure 8. The measured open charm 
Figure 7. The measured open charm spectrum compared with NLO pQCD 
spectrum compared with PYTHIA and the calculations for the charm quark. The 
Saturation Model. theoretical curves have been normalized to 

the measured cross section. 

the Saturation Model [19] and PYTHIA with parameters described in ref. [11] in Fig. (7). 
Note that the open charm spectrum has been scaled by Nb in = 7.5. One sees that 
the measured open charm spectrum is much harder than the PYTHIA prediction. 
However, the Saturation Model based on the kx factorization scheme does predict a 
harder spectrum in comparison with one obtained from the LO pQCD calculation in 
PYTHIA. 

In Fig. (8), the NLO pQCD predictions of the c-quark spectrum [20] are compared 
with the STAR open charm spectrum. The pQCD parameters in these calculations 
are m c (charm quark mass) = 1.2 GeV/c 2 , fip (factorization scale)=/iR (renormalization 
scale) =2 (transverse mass of the charm quark) for the MRST HO PDF curve, and 
m c =1.3 GeV/c 2 , fi F =fi R =m T for the GRV98 HO PDF curve. The theoretical curves 
have been normalized to the measured cross section though the cc cross sections from 
NLO calculations significantly underestimates the measured one (eg. by a factor of 
3 with the MRST HO PDF). It is interesting that the measured spectrum shape of 
the open charm coincides with those of the c-quark, especially the one obtained with 
GRV98 HO parton structure functions. Such a coincidence between the NLO c-quark 
spectrum and the measured open charm meson spectrum in charm hadroproduction 
has been observed in the fixed target experiments at low energies [21]. In order to fit 
with data, the intrinsic— kx model had to be introduced to counter-balance the effect 
of c-quark hadronization through a fragmentation function of Peterson form [1, 21]. 
At RHIC, however, adding intrinsic— kx with a moderate < k\ > would not help to 
overcome the momentum degradation of the fragmentation process because the charm 
quark spectrum is very broad. Such an observation may indicate that the hadronization 
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of a c-quark through the Peterson fragmentation function, which was established in 
charm photoproduction [22] and in e + e~ collisions [23], is not suitable in case of charm 
hadroproduction. Experimental data may suggest a fragmentation function more peaked 
towards the x — 1 region. A hadronization scheme of open charm production through 
quark recombination [24] was also proposed [25]. 
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